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ABSTRACT 

This thesis reports on the investigation carried out 
to study the effectiveness of sparse vector method in contin— ,, 
gency evaluation. The sparse vector method for solution of 
network equations has been examined in considerable detail, 
highlighting the scope for possible application. The sparse 
vector method has been tested for its viability in contingency 
evaluation based on Fast Decoupled Load Flow technique. 
Different sparsity oriented compensation methods have been 
reviewed and applied to contingency evaluation with a view 
to determining the suitability of sparse vector approach. 



CHAPTER 1 


INTRODUCTION 

1.1 GENERAL: 

Modern day power systems are all generally large and 
complex. Planning and operation of these systems require a 
number of studies to be carried out both in the off-line 
mode and on the real time basis. Most of these system studies 
necessitate solution of network equations involving large 
sparse matrices, which characterise the system. Sparse matrix 
methods [1,2] have been widely used for the solution of 
network equations, particularly from the view point of reducing 
computer storage requirement. However, for the system 

studies, to be carried out on real time basis, like contingency 
evaluation to ensure system security, it is important that 
the solution techniques be computationally efficient. In 
this context, solution techniques exploiting vector sparsity 
have been reported in the literature [3], both for planning 
and operation studies. 

1.2 SPARSE VECTORS: 

A system of linear algebraic equations may be 
represented as 
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where, 

[a] is the system matrix 
X and B are vectors 

In Eqn. (1.1), vector X is commonly referred to as the solution 
vector, while the right hand side vector B is called, the 
independent vector. In certain power system studies, either 
the independent vector B may be sparse in structure or the 
interest may be only in determining certain components of the 
solution vector X* In such a situation, vector B or vector X 
are termed as sparse vectors. 

In the commonly used method for the solution of linear 
algebraic equations, using factors of sparse A matrix, the 
network solution is obtained by performing full forward and 
full backward operations on the nonzero elements. However, 
with sparse vectors it is possible to avoid forward and 
backward operations on all the nonzero elements involved in 
the factors of sparse A matrix. The latter is termed as the 
sparse vector solution approach. 

In recent years* attempt has been made to use sparse 
vector solution approach for various power system studies. 
Brandwaj n and Tinney have applied the sparse vector 

approach to advantage in fault analysis [ The effectiveness 
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of this approach has also been investigated in case of 
sparsity oriented compensation methods by Alsac et.al. [5]* 
Possible application of sparse vector method in optimal 
power flow and economic despatch has been described in [6 ] 
and [ 7] respectively. 

Besides, certain other areas have also been identified 
for the application of sparse vector methods. These include 
contingency evaluation or outage studies, partial refactori— 
zation, parallel processing and network equivalencing [3 ], 

1.3 CONTINGENCY ANALYSIS: 

Contingency evaluation appears attractive from the 
point of view of application of sparse vector method . Conti- 
ngency evaluation which is a key function of security analysis 
is under taken on real time basis to analyse the outcome of 
contingencies such as line or transformer outages or loss of 
generation etc. on the security of the system. 

A number of fast-on line approximate algorithms [8,9^11] 
have been proposed for contingency evaluation. The techniques 
used fo.r contingency evaluation must possess the essential 
features of fast computational speed, reliability, reasonable 
accuracy and moderate memory storage requirements. The 
present trend in contingency evaluation algorithms is to 
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compensate the base case load flow results, using certain 
matrix modifications which simulate the effect of outages. 

This approach is preferred compared to the refactorization 
approach in which the base case system matrix is modified 
and then refactorized to obtain the network solution for the 
post-contingency operating state [5 ], 

It is essential to highlight that in contingency 
evaluation, matrix modifications to reflect the configuration 
changes arising due to outages are derived on the basis of 
well known Inverse Matrix Modification Lemma [5, 12, 12] .Since these 
methods of matrix modifications, offer a basis of providing 
compensation for network changes, these are generally 
referred to as compensation methods. A wide range of compen- 
sation methods has been reported in the literature [2,5], It 
is in these compensation methods for contingency evaluation 
that sparse vector method may find possible application. 

1.4 OBJECTIVE AND SCOPE OF THE THESIS: 

From the review of available literature, it seems 
important to further explore the potential of sparse vector 
approach as regards its application, merits and limitations. 

The objectives of this thesis, thus, are: 
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a) To study the characteristics of sparse vectors and 
to examine the sparse vector method in detail with 
particular reference to its implementation fox solu- 
tion of linear algebraic equations. 

b) To investigate the possible application of sparse 
vector method for contingency evaluation employing 
compensation techniques. 

1*5 OUTLINE OF THE THESIS*. 

Chapter 2 deals with the study of sparse vectors as 
regards their characteristics. The derivation of sparse 
vector method for solution of network equations and its 
implementing algorithms are described. Further, the effi- 
ciency of sparse vector approach is examined in obtaining 
solutions of linear algebraic equations by conducting 
several, case studies. 

The application of sparse vector solution method for 
contingency evaluation is described in Chapter 3, The 
contingency evaluation uses compensation technique, to 
simulate the effects of outages. The effectiveness of the 
sparse vector approach is investigated through case studies 
on a number of sample power systems. 
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With reference to contingency analysis a comparative 
evaluation of various compensation techniques, considering 
the application of sparse vector method, is reported in 
Chapter 4* The computational requirements of various 
compensation techniques are compared through case studies 
on different sample power systems. 

Chapter 5 reviews the major contribution of the thesis 
and suggests the scope fox future work. 



CHAPTER 2 


SPARSE VECTOR METHOD 

2.1 GENERAL: 

. While the sparse matrix methods for the solution of 
large scale power system problems are well established, vector 
sparsity has not been fully exploited despite its advantages. 
Tinney et.al. [3 ] have highlighted the usefulness and 
potential of the sparse vector approach for the solution of 
linear algebraic equations. An attempt has, therefore, been 
made in this chapter, to study the effectiveness of sparse 
vector methods in the solution of linear algebraic equations. 

2*2 SPARSE VECTOR APPROACH: 

Consider a system of linear algebraic equations 

[A] X = 8 (2.1) 

where 

[A] is -the system matrix of order N 
X is the vector of variables or solution vector 

B is the independent vector 

Sparse vector methods require the system matrix [A] to be 
factored as 


[A] = [L] [D] [U] 


( 2 . 2 ) 
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where, 

[Lj is the lower triangular matrix with unity diagonal 
[D] is the diagonal matrix and 

[U] is the upper triangular matrix with unity diagonal. 

However, sparse vector methods are equally applicable to the 
cases where [A] is factored as 

[A] = [L] [U] (2*3) 

The matrix U is the transpose of matrix L when [a], is 
symmetric, which is normally the case in respect of power 
system admittance matrices. 

The solution vector X using triangular factorization 
can now be obtained from Eqn. (2.1) and Eqn* (2.2) as 

[L] [D] [U] X = B (2.4) 

Substituting [U] X = Y and [L] [D] = [L* ], equation 
(2.4) gives 

[L l ] Y = B (2.5) 

or Y =* [L* J ^ B (2.6) 

Solution of Eqn. (2.5) gives an intermediate vector Y, which 
in turn is used to determine the solution vector X from 

[U] X = Y (2.7) 
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Then, 

X = [U]" 1 Y (2.8) 

Equation (2.6) defines a sequence of forward substi- 
tution operations on [L* ], while Eqn. (2.8) defines a 
sequence of backward substitution operations on [u]. The 
operations can be performed by rows or columns. An analysis 
of the algebra involved in the sparse vector approach will 
show that for sparse vector method , it is essential that the 
forward solution of equation (2.6) be performed by columns 
and in the backward solution of Eqn. (2.8), the elements of 
[U] be accessed row-wise [3 ]. 

2*2.1 Fast Forward/Fast Backward Operations: 

When the independent vector B is full, the sequence of 
forward substitution operations on elements of [ L* ], to obtain 
intermediate vector Y is required to be carried out for each 
column, in the ascending order of columns. This may be 
termed as the 'Full' Forward Method*. If the vector B is 
sparse, only a subset of the columns of [L* ] is needed for 
the forward solution of Eqn. (2.6). This is called the 'Fast 
Forward* (FF) solution* 

Likewise, if all the elements of the solution vector 
X are required to be computed, it calls for 'Full Back' 
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substitution, in which case elements of each row of [u] 
are accessed* However, when the interest lies in knowing 
only certain elements of the vector X* then only a subset of 
rows of U is needed for the solution of Eqn. (2.8). This 
is termed as the 'Fast Backward’ (FB) solution. 

For the application of sparse vector method , it is 
necessary to efficiently identify the subsets of columns of 
[L’J and [U] for FF and FB solutions. The subset of columns 
for FF solution is a function of the sparsity structures of 
[L’ ] and B. Similarly, the subset of rows for FB solution 
is a function of the sparsity structures of [U] and X. The 
sparsity structures of [L* ] and [U],in turn, are the functions 
of the sparsity structure of system matrix A. 

Some algorithms have been proposed to identify the 
subset of columns of [L’ ] for Fast Forward solution and the 
subset of rows of [u] for Fast Backward solution [3 ]. 

These algorithms employ the concept of 'factorization paths’, 
which enables identification of the relevant subset of 
columns and rows. 

2.2.2 Factorization Paths: 

Before defining a factorization path, it would be 
appropriate to describe a ' singleton’ , because the path 
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concept has been developed over the ' pa th associated with a 
singleton [3], A vector with only one nonzero element 
may be referred to as a singleton* A* factorization path* 
or ’path’ can be defined as an ordered list of columns of 
[L‘ ] for Fast Forward operations or rows of [u] for Fast 
Backward operations. The path for a singleton may be deter- 
mined from the following algorithm. 

a) Search for the nonzero element in the vector B and 
let k denote its location in the vector. Record k 
as the first element in the path. 

b) Move to the column k of [L’ ] or row k of [Uj and 
determine the location of the first nonzero element 
in column k of [L* .3 or row k of [uj. In searching 
for the first nonzero element, the diagonal elements 
of [L* ] and [U] are not to be considered. 

c) Replace the previous value of k in (a) with the 
location of the first nonzero element traced in (b) 
and list it in the path. 

d) Repeat the steps (b) and (c) till k = N (N being 
the order of [L*] or [u3). 

It needs to be mentioned that a path is executed in forward 
order for FF solution and in reverse order for FB solution. 
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The formation of the path associated with a singleton 
is illustrated through an example. Consider a 10 node 
network for which the structure of [£] is shown in Fig. 2.1* 
While the star ('*) signs in the figure indicate the locations 
of the nonzero elements initially present in [a] , the 00 signs 
represent the location of additional nonzero fill-in elements 
created in obtaining the factor [L* ]. Now, for a singleton 
B with nonzero in location 1, using the preceding algorithm, 


the path is obtained 

as 
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When matrix A is symmetric, its factors [L* ] and [u] 
are also symmetric. In such a case, the path associated with 
a particular element in the vector X for FB solution would be 
the same as the path determined for FF solution, for the same 
location of nonzero element in vector B, Nevertheless, the 
subset of rows of [u] or path for FB substitution can be 
determined in an identical manner, using the preceding 
algorithm, irrespective of whether [L* ] and [U] are symmetric 
or not. The sparsity structure for a non- symmetric matrix 
U is shown in Fig, 2.2, For a singleton with nonzero at 
location 2 i.e. in case the interest lies in obtaining 
only the second component of vector X, the path may be traced 
as 2,3,6,8,9,10 . 

However, as mentioned earlier, it is essential .that 
the path for Fast Backward solution be executed in reverse 
order implying that in case ,the second component of vector 
X is required to be determined, then to use the Fast Backward 
solution, the path has to be executed in reverse order viz. 

10,9,8,6,3,2 • If matrix A is stored exploiting sparsity 

through indexing arrays, the path for a singleton can 
also be determined directly from the indexing arrays, 
instead of searching and testing for the nonzero entry, as 
carried out in the above example, where the storage of full 
A matrix is assumed. 



14 


NODE 

1 

2 3 4 5 6 

7 8 

9 

10 

1 

* 


* 

* 

* 

2 


* * 




3 


* * 




4 


* 

oo 



5 


* 

* 


* 

6 


* 

oo 

* 


7 






8 



* 

* 


9 




* 

OO 

10 





* 

Fig . 

2.2: 

Sparsity structure 

of matrix 

u 



In a general situation, it is possible that the vector 
B may have more than one nonzero elements or the interest 
may be in evaluating more than one component of the vector X, 
Corresponding to each nonzero entry in the vector B, a path 
can be obtained as explained earlier in case of a singleton. 
The union of these individual paths defines the complete 
factorization path for the sparse vector consisting of more 
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than one nonzero entry* The algorithm for obtaining the 
complete factorization path for a general sparse vector is 
described below [3], 

i) Search for the first nonzero element in the vector B 
and trace its singleton path as explained earlier. 

This path segment should always end with node N. 

ii) Identify the next nonzero node, which is not already 
in the partially completed path and determine its path 
in the same manner, until a node already in the partially 
completed path is encountered. 

iii) Next, join this newly formed path segment to the 
partially completed path list. 

iv) Steps (ii) and (iii) are repeated till all nodes in 
the sparse vector are in the path, 

v) The path list is arranged either in ascending order 
for Fast Forward solution requirement or in reverse 
order for Fast Backward solution. 

Based on the above algorithm, the determination of 
the complete factorization path for the FF solution is 
described for the sample 10 node system example. Consider a 
vector B with nonzero elements at location 1,3,6 and 7. 

Using step (i), the path corresponding to nonzero in location 1 
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for the sparsity structure of matrix L* represented in 
Fig, 2.1 is traced as 1,4,5,8,9,10 . The next nonzero being 
in location 3, the second path segment using step (ii) is 
obtained as 3,7 . Step (iii) gives the partially completed 
path as 3,7,1,4,5,8,9,10 . For the nonzero in location 6, 
the path segment is determined as 6 , comprising only one 

node. The partially completed path is then enlarged to give 
6,3,7,1,4,5,8,9,10 • Path segment corresponding to nonzero 
in location 7 is not required to be traced, as node 7 already 
appears in the computed path list. The elements of the path 
list are then arranged in ascending order as 1,3,4,5,6,7,8,9,10 
to give the factorization path for the sparse vector B for 
Fast Forward solution. 

A pictorial representation for the sparse vector with 
4 nonzeros as discussed is shown in Fig. 2.3. The path for 
node 1 is formed by tracing a path from this node in ascending 
node sequence, to the last node (node 10 in this case). The 
partial paths emanating from nodes 3 and 6 are traced out 
similarly, until they meet the already formulated path at a 
junction (node 8). 

No additional path is required to be generated for 
node 7, as it is already contained in the path originating 
from node 3. The graph of Fig. 2.3 containing the union of 
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these path segments describes the subset of columns of [L* ] 
for FF solution. An analogous approach would be valid for 
obtaining the subset of rows of [U] for FB solution. 

2.3 TEST EXAMPLES: 

To illustrate the sparse vector approach described 
in the previous section, the solution of a system of linear 
algebraic equations corresponding to 20 node, 40 node and 
80 node synthetic network is obtained employing 

a) Fast Forward and Fast Backward solution 

b) Full Forward and Full Backward solution 

Computer programs have been developed for both the 
solution methods. In the programs developed, full storage 
of A matrix is contemplated and as such optimal ordering 
is not attempted. The effectiveness of the sparse vector 
solution approach for the various sample systems is investi- 
gated through 

i) variation in the number and location of nonzero 
elements in the independent vector B. 

ii) Variation in the degree of sparsity of the 
system matrix A. 
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With the variation in the number of nonzero elements 
in the vector B, the computation time requirement for the 
solution of a system of linear algebraic equations corres- 
ponding to a 20 node, 40 node and 80 node network is tabulated 
in Table 2.1. The computation time reported includes the 
time for the factorization of matrix A. 

It can be seen from Table 2.1 that the computation 
time for FF/FB solution is largely dependent on the number 
of nonzeros in the B vector. This fact is particularly 
obvious in the solution of the 80 node system. Furthermore, 
for the three cases studied, with large number of nonzero 
entries in the B vector, the Full Forward/Full Backward 
solution is more beneficial. This is because of the fact 
that as the sparsity of vector B decreases i.e. as the number 
of nonzero elements in vector B increases, the path finding 
algorithm becomes increasingly inefficient, thereby accoun- 
ting for the increased computing time of FF/FB solution 
method. 

However, when the B vector is reduced to a singleton, 
the solution is obtained faster with the FF/FB approach. 

This observation too is very evident from Table 2.1, in all 
the three cases. Moreover, for a singleton B, reduction in 
execution time for FF/FB solution is quite appreciable 



Table 2.1: Computational Timings for Solution of 20 Node, 40 Node, and 
80 Node System Equations, 
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as the size of the system is increased. For a smaller 
system, the Full Forward/Full Backward solution is more 
relevant. It may be worth mentioning that case studies 
reported in [3,5]., to illustrate the effectiveness of sparse 
vector method have been carried out on systems having 
400 to 2000 nodes. 

During the course of the study in Table 2.1 it was 
also noticed that the variation in the location of the nonzerc 
element, in a singleton B did not have any significant effect 
on the computation time, A similar observation was also made 
when the magnitude of nonzero element was changed, keeping 
its location fixed. 

The effect of variation of degree of sparsity of the 
matrix A, on the execution time for 40 node and 80 node 
system is shown in Table 2,2, The sparsity of matrix A has 
been reduced to in both these cases, 

A comparison of Table 2,2 with Table 2.1 for the 40 
node system reveals that with the reduction in the sparsity 
of [A], the time required for FF/FB solution reduces drasti- 
cally, while there is no appreciable change for the time 
required for the Full Forward/Full Backward solution, even 
with large number of nonzeros in vector B. The reduction 



Table 2.2: Effect 'of Variation of Degree of Sparsity on 40 Node and 
80 Node Systems. 
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in the nonzero elements in the A matrix leads to reduction 
in the path length, consequently resulting in faster 
computing time, A similar trend is also observed when the 
80 node systems of Tables 2.1 and 2.2 are compared. In 
general, in both cases, | FF/FB solution may be preferred when 
vector B is highly sparse. 

Computer programs developed for the solution using 
FF/FB approach are also capable of solving a system of 
equations for different independent vectors B. The program 
avoids the need for refactorization of matrix A, with every 
Change in B vector. 

Table 2.3: Comparison of FF/FB Approach and Full Forward/ 
Full Backward Approach for a 80 Node System. 

(System Matrix £ has 
byi nonzero elements). 


SI. 

No. 

No. of 
nonzeros 

Execution Time 
( in msecs) 

FF/FB 

solution 

Full Forward/Full 
Backward solution 
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80 
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141 

2 

56 

241 


141 

3 

2 

116 


146 

4 

1 

105 


169 

5 

1 

98 


141 
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The time requirement for both the solution approaches 
is recorded in Table 2.3. The execution time shown does not 
include the time for factorization of A. Once again, it 
becomes obvious from the results of Table 2.3, that with a 
small number of nonzero elements in vector B, FF/FB solution 
method is preferable. 

2.4 CONCLUSION: 

This chapter reports on the use of sparse vector 
approach for the solution of a system of linear algebraic 
equations. The effectiveness of the FF/FB solution method 
employed in the sparse vector approach is compared to the 
normal Full Forward/Full Backward solution method, through 
various case studies. It is observed that the sparse vector 
methods can be preferred in situations where the system 
under study is extremely large and the independent vector 
is highly sparse, if not a singleton. 



CHAPTER 3 


CONTINGENCY EVALUATION USING SPARSE VECTOR METHOD. 1 

3.1 GENERAL: 

With the growing complexities of the present day- 

power systems, security analysis has become an important 

aspect. One of the functions of security analysis which is 

commonly known as contingency evaluation is to determine 

whether the normal system is secure or insecure. Based on 

this information, the preventive action to be taken is 

determined. In order to ensure satisfactory system operation, 

* 

contingency evaluation is generally carried out in real time. 
In view of this, it is important that the techniques for 
contigency evaluation should be computationally efficient. 

In this context, this chapter describes the application of 
sparse vector methods to contingency evaluation. 

3.2 CONTINGENCY EVALUATION: 

A power, system generally undergoes contingencies 
like outage of transmission line, outage of a transformer 
or loss of generation etc. Investigation of the effects of 
these outages is essential both during system planning and 
operation. Outages may be grouped into two categories. 
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a) Network outages 

b) Power outages 

A network outage involves changes in the network 
admittance parameters and includes the following 

i) outage of a transmission line 

ii) outage of a transformer 

iii) outage of a shunt capacitor or reactor 

A power outage brings about changes in bus power injections 
and includes the following 

i) Loss of generating unit 

ii) Sudden loss of a load 

iii) Sudden change in flow in an. inter tie. 

It is not practicable to test a power system for every 
conceivable contingency.. Based on historical records and 
operating experience,, a nominal contingency list is prepared* 
These contingencies are then processed in an order that 
takes into account their probability and severity* However, 
it is has been observed that of all the contingencies, the 
transmission line outages are quite common and thus attention 
is confined in this chapter to network contingencies. 

Network outages result in a change of system configu- 
ration. -One of the obvious ways for determining the effect 
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of these contingencies is to carry out fresh load flow for 
the changed system configuration. But carrying out complete 
load flow for every contingency is both time consuming 
and costly, especially if many cases are to be studied. The 
alternative is to simulate the effect of these outages to 
obtain corrections, which may be used to update, the base 
case load flow solution, to represent the post contin- 
gency operating state of the system. 

Several methods for contingency evaluation based on 
linearized a.c. load flow, distribution factors, d.c. load 
flow,, full a.c. load flow etc. are being used *by various 
pow'er utilities [8-11,13-15], Fast Decoupled Load Flow [8] 
provides a straightforward method for contingency evaluation. - 
An iteration cycle of Fast Decoupled load flow consists of a 
phase angle iteration followed by a voltage magnitude iteration. 
After obtaining the base case load flow for the normal 
system, each line outage case in the list of contingencies is 
then processed by the compensation method using a *9* itera- 
tion to detect probable overloads. If voltages are a 
problem, then 1 9* iteration is followed by a 'V' iteration. 

Based on Fast Decoupled load flow method, the base 
case load flow is obtained using the following equations. 
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[ A©] and [A.V 1 ] are solution vectors for phase angle 
corrections and bus voltages corrections. 

Matrices [B f ] and [B 11 ] are derived from Y Bus matrix. 

Both [B«] and [B * * ] are real and sparse. Since they contain 
only network admittances, they are constant and need to be 
evaluated once only at the beginning of the study. 

Network outages like line or transformer outages can 
be simulated using the Inverse Matrix Modification Lemma [5 ]. 
Outages are reflected in base case matrices B’ and B' * by 
modifying the affected elements in row k and row m. For 
each line or transformer outage, only four elements of the 
matrix B f and B ’ 1 undergo change. Without actually modi- 
fying B' and B ’ 1 and then performing a repeat load flow, a 
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row vector M which simulates the effect of disconnected 
line or transformer is introduced. In case, the line 'km 1 
between the kth and mth nodes gets disconnected , the row 
vector M has the following structure. 

kth 51 th 

co lumn co lumn 

M = [000 a 0 00-1 000 ] 

(3.3) 

whe re , 

a = off-nominal turns ratio of the transformer 
between buses k and m 
= 1 , for a line 

The matrix Bicorresponding to the post-contingency 
configuration of the system, consequent to the outage, is 
obtained as 

[Bj_] = [B‘] + b. [M T ][M] (3.4) 

where, 

b = line or transformer series admittance between 
nodes k and m 

Then, 

[ B ,*]" 1 = |[ B ' ] + b [M T ] [M] f " X (3.5) 

Using the Sherman-Morrison formula [12] and expanding 
Eqn. (3.5), we have 
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[B| r 1 = [B ' ]"" 1 -b j l+b[M][B’]~ 1 [M r ]f ]~ 1 [M T ] [mJ [B 1 ] -1 j 

° r , 

[B|1=[B* ]“ 1 ~b 1 b(^ +Cm][B*]- 1 CM T ]]}~ 1 [B»]“ 1 [M T ][M][B f 3’ 1 
or, 

[B*]= [B ’ ] _1 -He +[ M ][B , ]“ 1 [M T 3)i'" 1 {[B*]” 1 [M T ][M][B»]“ 1 | 

x l u o L 3 

(3.6) 

Let X = [B> ] _1 M T (3.7) 

Substituting Eqn. (3.7) in Eqn, (3,6) 

[B']” 1 = [6*]" 1 - j A + M X) “ 1 X M[B’ j” 1 (3.8) 

1 £ D “) 

Replacing + M X [■ = C, Eqn. (3.8) becomes 

/ U ) 

[B^r 1 =[B«]“ 1 - c X M [B']” 1 (3.9) 

Phase angle corrections for outages axe now derived 
as under. 

Equation (2.1) for the post-contingency state 

yields 

[B-^3 [ A ©}_] = [ -j^y] (3.10) 

or, 

[4©^ = [B^ -1 [ “ |Y[] (3.11) 
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Substituting for from Eqn. (3,9) in 

Eqn. (3,11) 

[ -Ae^ ^KbT 1 - c X M [B']- 1 ! C-pf- ] 

or, 

[ ae 1 ] = [B *]’ 1 cx m [B-rHf^-] 

But [B.] [^P_] = [a 9 ] 

Hence [ A © x ] = [ A© ]- CXMUe] 

Equation (3,14) gives the correction vector as 

00 = -c x m[aq] 

In exactly similar fashion, bus voltages corrections for the 
post outage configuration can be computed. 

After each solution of Eqn. (3.1), [A 9] is corrected 
by an amount 

S'Q = -C X M [ A ©] 

Similarly, after each solution of Eqn. (3,2) [AivQ is corrected 
by an amount 

6 V = -C* * X' ' M' 1 |>V] (3.16) 

where superscript * ’ is used to denote vectors and variables 
used for obtaining bus voltage corrections. 


(3.12) 

(3.13) 

(3.14) 

(3.15) 
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The technique outlined above, based on Fast Decoupled 
load flow method can be used for multiple outages as well. 

For multiple simultaneous outages, the above procedure is 
applied recursively. The base case solution vector is 
corrected successively, as the effect of each branch outage 
is introduced one at a time. For a set of 1 n* outages, the 
recursion is 


[ Ae i ] = [ Wit'-i-d 

for i = 1,2, .. .n 


(3.17) 


The scalars C. and the vectors X. are calculated at the 

X X 

beginning of the case from the recursive algorithm 


x i =[BT 1 «i T 


(3.18) 

X i = x ± -c for j = 

i— 1 ) 

(3.19) 

C j _ = (l/b j _ - M j _ X.)- 1 for 

3 1 1 i= 1*2.* • • 

(3.20) 

■ recursive procedure is a 

pplied, when bus 

voltage 


corrections [ t iVfjare to be computed, for multiple simul- 
taneous outages. 

It is of interest to point out that vector M in 
Eqn. (3.7) is highly sparse having only two nonzero elements. 
The sparsity of this vector can be exploited using sparse 
vector method in the solution of Eqn. (3.7). 
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3.3 ALGORITHM FOR COOT INGE MCY EVALUATION: 

The post contingency state of the system using Fast 
Decoupled load flow technique is obtained using the algorithm 
described below. The algorithm proceeds with the assumption 
that the base, case load flow solution is available. 

a) B* matrix is obtained from Y Bus matrix in the 
normal manner. 

b) Mismatch vector xP/jV| is computed 

c) Row vector M is formed to simulate outages and 

T T 

M is then obtained. This vector M is highly 

sparse. 

d) Vector X is computed using Eqn. (3.7). In obtaining 
vector X, triangular factors of B* are employed, which 
are computed once i n the beginning. 

e) Correction vector £'&. is next computed from Eqn. (3.15)., 

f) Updated A'G^is then obtained from ^qn. (3.14) and 
consequently bus angles corresponding to post-contin- 
gency state are computed. 

g) Analogous steps as in (a) to (f) are also performed 
for computing updated bus voltages pertaining to the 
post-contingency state of the system. 

Generally one iteration of * 9 ' and one iteration of 'V' axe 
performed, as this gives reasonable idea about the system 
security Cl3]« 
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3.4 CASE STUDIES: 

The application of sparse vector methods in contingency- 
evaluation is illustrated through outage studies on the 
following systems. 

a) IEEE 14 Bus system 

b) IEEE 57 Bus system 

c) Practical 100 Bus system of a Utility 

Computer programs based on Fast Decoupled load flow 
have been developed for contingency evaluation employing (a) 
Fast Forward/Fast Backward approach and (b) Full Forward/ 
Full Backward approach for the solution of the system 
equations. 

The computation time for various segments of the 
algorithm using the two approaches is tabulated in Table 3.1 
for one iteration of phase angle correction and one iteration 
of bus voltage correction. 

From Table 3.1, it Is evident in all the cases, that 
while Full Forward/Full Backward approach is faster for 
computation of uncorrected £v0 and uncorrected A jVi, the 
FF/FB approach is beneficial for the solution of vectors X 
and X’ ’.On the basis of this observation, a hybrid approach 



Table 3.1: Computational Performance of FF/FB Approach and Full Forward/Full 
Backward Approach in Contingency EvaUftation. 
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involving Full Forward/Full Backward approach for obtaining 
uncorxected a© and uncorrected and FF/FB approach for 
computing vectors X* and X' * was also examined. The time 
required for computation in this case is recorded in Table 3,2, 

Table 3,2: Analysis of Hybrid Approach in Contingency 
Evaluation 


SI, Program 

No* segment 

Computation 
( in msecs) 

Time 


IEEE 14 Bus system 

i£e£ 57 Bus 
system 


1 Uncorrected [a© ] 

17 

854 

4480 

2 Uncorrected [a IV 1 ] 

0 

586 

2499 

3 X* 

6 

801 

4366 

XI 

9 

562 

2385 

J5 Complete 

1 

285 

5181 

22627 


Comparison of computation timings in Table 3,1 and 
Table 3,2 indicates that it would be beneficial to use 
hybrid approach instead of Full Forward/Full Backward approach 
or FF/FB approach alone, for contingency evaluation. It is 
expected that the efficiency of the hybrid approach would be 
predominantly better than the other approaches, with increase 
in the system size. 
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3.5 CONCLUSION: 

Contingency evaluation is now one of the very important 
studies in the planning and operation of a power system. It 
forms a desirable part of a comprehensive system security 
monitoring process. As commercial power system utilities 
these days, evaluate contingency performance of the system 
very frequently, computationally efficient methods for conti- 
ngency evolution are greatly sought after. In this chapter, 
fast decoupled load flow technique exploiting vector sparsity 
has been examined. It is observed that a hybrid approach 
involving a judicious mix of FF/FB approach and Full Forward/ 
Full Backward approach would improve the computational 
efficiency. 



CHAPTER 4 

COMPARATIVE STUDY: ' OF COMPENSATION TECHNIQUES FOR ’ 

CONTINGENCY EVALUATION 

4.1 GENERAL: 

Compensation methods can be used advantageously for 
almost any application in which each new case requires a 
small modification in the network matrix with respect to the 
base case. Several compensation methods have been reported 
in the literature [5 ], These compensation methods d.iffer 
from each other in the manner in which compensation is applied 
This chapter deals with application of sparse vector 
approach for contingency evaluation based on different 
compensation techniques, 

4.2 COMPENSATION METHODS: 

Compensation methods, which are based on Inverse 
Matrix Modification Lemma are used to simulate outages and 
obtain post-contingency state of a system, subjected to 
outages, A brief outline of the compensation methods is 
discussed in the following. 

In the base case load flow, Fast Decoupled load flow 
equation for obtaining phase angle correction vector A ©, 
as already stated in preceding chapter is 
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[B- ] [ AS] = [-^r 3 (4.1) 

where symbols have their usual meaning* 

For changes in the system configuration due to outages, 
Eqn. (4,1) gets modified to 

[ B l] [Zi0 1 J = [-|£j-] (4*2) 

The solution of Equation (4.2) gives the post- 
contingency phase angle as 

[ A e x 3 = [B 1 ]" 1 [ -|y|" 3' (4,3) 

Applying Inverse Matrix Modification Lemma to Eqn, 

(4,3) and proceeding in the same manner as for Eqn. (3.5) 

[ 40 1 ]= |[B ! ]~ ;L -C[B»]“ 1 M T .M[B» J” 1 ! [“J7J-3 '(4.4) 

where again symbols used have their usual meaning* 

Substituting M T 

[ 4.0 1 ] =j[B»r 1 -C X.M [B* r 1 f (4.5) 

Equation (4,4) can be arranged to give different computational 
arrangements, for determining phase angle correction vector 
Taking [B 1 over to the right in Eqn. (4.4) 

[ *6^ =<I-C[B' ] -1 M T M][B'] _1 (4.6a) 

where I is identity matrix. 

Alternatively, moving [B 1 ] -1 to the left in Eqn. (4.4), 
vector is obtained as 

[ A© x ] = [B* r 1 | I-C M T M[B> r 1 ! [-JVf] ( 4 * 6b ) 

Another arrangement is derived when [B*] - ^ is replaced by 
its factor [U]" 1 [Lj"’ 1 , 


so as to give 
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[ ^3 = [u]"* 1 Ii-Cl’T 1 c mVEuF 1 ! Cl]" 1 C-jyp] 

(4.6c) 

The different computational arrangements of Eqn. (4.6a), 
(4.6b) and (4.6c) provide the basis on which different compen- 
sation methods have been evolved. Each of these equations can 
be analysed to have a central calculating step in which 
network solution is performed by forward and backward substi- 
tution, using the sparse [L] and [u] factors of matrix B*. 
Additionally, there is a compensating step in which compen- 
sation for network modification is applied to the system. 


Compensation methods, are classified depending on the 
stage at which compensation for network modification is 
applied to the system. Vi/hen compensation is provided for, 
after the network solution is carried out initially^as per 
Eqn. (4.6a), the method is referred to as post-compensation 

/"\ 15 

method. In Pre-compensation, the independent vector W 
is compensated before the network solution as per arrangement 
of Eqn. (4.6b). Equation (4.6c) gives the equation for Mid- 
compensation, where the intermediate vector in the middle of 
the 'network solution i.e. after the forward solution but 
before the backward solution, is compensated. 


The basic steps involved in the three different 
compensation methods, for computation of phase angle 
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corrections are outlined in Table 4.1. Analogous steps would 
be applicable for determining bus voltage corrections. 

It will be seen from Table 4*1 that column vector 
which is highly sparse is being used for computation of the 
following vectors,. 

i) X, for Post-compensation from equation (4.8) 

ii) Xpfor Pre-compensation from equation (4.13) 

iii) W and W^»for Mid-compensation from equations (4.20) 
and (4.21) respectively. 

This provides the scope for examining the application of 
sparse vector approach in the compensation methods discussed. 
In fact, it may be noted that the method used for contingency 
evaluation in the previous chapter relies on the post-compen- 
sation technique. This can be easily inferred by comparing 
Eqns. (4,5) and (4.6a) with Eqn. (3,12). 

4.3 CASE STUDIES: 

Case studies reported in the previous chapter dealt 
with application of sparse vector methods in contingency 
evaluation based on post-compensation technique. Sparse 
vector methods may be applied to other compensation methods 
as well, because of the obvious sparsity structure of M and 
M 1 vectors. Separate computer programs incorporating 
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TABLE 4.1 

BASIC SOLUTION STEPS FOR POST, PRE AND MID-COMPENSATION 


POST-COMPENSATION 

. _ ; 

SI. No. 

STEPS 

1 

Obtain network solution 



4e = [b- r 1 [^-] 

(4.7) 

2, 

Compute compensation vector 



X = [ B * 3 ~ 1 M T 

(4.8) 


Z = M X 

(4.9) 


; C = (1/b + Z]"" 1 

(4.10) 


00 = -C X M A© 

(4.11) 

3. 

Apply compensation 



A 0 1 = 4© + £e 

(4.12) 


Table 4.1 contd. 
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PRE-COMPENSATION 


SI. No. 

STEPS 

i 

Compute compensating vector ( OP) 



X, = [B' 1 ]- 1 M T 

(4.13) 


h =*i « T 

(4.14) 


Cjl = (1/b + Z)" 1 

(4.15) 


[ 5p ] - - c i-m T Xi T 

(4.16) 

2. 

Apply compensation 



compensated [p^j-] = [p^p] + £ '■ ^ 3 

(4.17) 

3 

■i 

Obtain network solution 


• ] 

i 

"i 

£6-^ = [B'T 1 j compensated [p^|*] \ 

(4.18) 


Table 4,1 contd. 
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iViID—COJvlPE NSAT I ON 

SI. No, 

STEPS 

1 

Start network solution 



F a = [L] -1 [ 

(4.19) 

2. 

Compute compensating vector AF 



w = [l]" 1 • M T 

(4.20) 


w x = (U 1 )" 1 . m t 

(4.21) 


z = W X T . w 

(4.22) 


C = (1/b + Z)" 1 

( 4 * 23 ) 


F = -C.W.WjJ F x 

(4.24) 

3 

Apply compensation 



F = ?i + F 

(4.25) 

4 

Complete network solution 


i 

Ae 1 = [U]" 1 . F 

(4.26) 
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Fast Forward/Fast Backward (FF/FB) approach and Full Forward/ 
Pull Backward approach were developed to examine the relative 
efficiency of sparse vector methods in the following compen- 
sation techniques. 

a) Pre-compensation 

b) Mid-compensation 

Pre-compensation and mid-compensation have been tested 

in contingency evaluation, which is based on Fast Decoupled 
load flow method. 

-The computation time required for various program 
segments and the complete program, for Pre-compensation and 
mid— compensation techniques is shown in Tables 4,2 and 4.3 
respectively. In Table 4.2, computation time for Pre-compen- 
sati.on technique has been recorded for the following 
pr-ogram segments, 

i) Phase angle correction from Eqn. (4, 18), where 

the independent vector is non-spar se. 

ii) Bus voltage correction using equation analogous 

to Eqn. (4.18) . 

iii) Vector fusing Eqn. (4.13), where the independent 
vector M T is highly sparse. 

iv) Vector using equation analogous to Eqn. (4.13). 



Table 4.2: Computational Performance of Pre-Compensation Techniqu 
Contingency Evaluation. ' 
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Also, the program segments computed in Table 4.3, for 
Mid-compensation technique are 

i) Vectors F^ and A©^ for phase angle corrections, 

using equation (4.19) and (4.26) respectively, where 
the independent vectors are non-spar se, 

ii) Vector F^y and AVj^ for bus voltage corrections, 

using equations analogous to Eqns. (4.19) and (4.26). 

iii) Vectors W and from equations (4.20) and (4.21) 

T 

respectively, where column vector M , is highly 
sparse. 

iv) Vectors Wy and W^y from equations analogous to 
Eqns. (4.20) and (4.2l), 

A study of the computation timings for various program 
segments in Table 4.2, for Pre-compensation shows that 
neither of the two approaches i.e. the FF/FB or the Full 
Forward/Full Backward approach is beneficial by itself. For 
all the sys tems studied, FF/FB approach exploiting vector 
sparsity is advantageous in the solution of equations for 
highly sparse independent vectors, whereas the Full Forward/ 
Full Backward approach may be preferred, when the independent 
vector is full or non— sparse. It will be seen that this is 
in concurrence with the observations made in the preceding' 
chapter, in the context of post-compensation based contin- 
gency evaluation. 
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An identical observation can be made from the analysis 
of computation time requirement recorded in Table 4.3, for 
contingency evaluation based on Mid-compensation technique. 
Here again, FF/FB approach is suited for highly sparse inde- 
pendent vectors, while Full Forward/Full Backward approach may 
be pursued when the independent vector is not highly 
sparse. 

In view of the above observations, again a hybrid 
approach combining both the approaches i.e. the FF/FB and the 
Full Forward/Full Backward approach,, to exploit their indivi- 
dual advantages, has been explored. The hybrid approach has, 
been applied for contingency evaluation, based on the 
three compensation methods under study and the computation 
timings are reported in Table 4.4. 

A comparison of Table 4.4 with Tables 3.1, 4.2 and 4.3, 
for each of the compensation techniques, brings out that the 
hybrid approach is more beneficial than either of the other 
two approaches, when the latter are used in isolation. 

It is also possible to make a comparison of the three- 
compensation techniques as regards the overall computing 
time from the results of Table 4.4. For the specific case 
of one outage at a time, the results in Table 4.4 indicate 



Table 4,4! Comparative Evaluation of Hybrid Approach in Post,pre and Mid-Compensation 
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8 W 1V 2354 

Complete; 22627 Complete 23287 9 Complete 43460 
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that post-compensation seems to be least time-consuming, 
though Pre-compensation and Post-compensation do not differ 
significantly in terms of computation time* 

However in the cases studied, Mid-compensation takes 
more time than the other two compensation techniques, though 
in reference [5 ], it has been reported that for a general 
power network problem, Mid-compensation is usually more 
efficient and faster as compared to the other two compensation 
approaches. 

A direct comparison of the computing time of the 
compensation methods given in Tables 4.2, 4.3, 4.4 and that 
reported in [5 ] is not appropriate. In reference [5 ] the 
formulation, has been divided into two phases. 

a) The Preparatory Phase - Here the calculation of 
certain terms of the inverse of matrix B* and B’*or 
their factors, to obtain Z and C is carried out. 

b) The Solution Phase - This comprises a standard 
network solution using the factors of the original un- 
modified matrices B' and B* 1 , and the compensating 
steps that introduce the effect of the modification. 

In the present case, no distinction has been made in 
the two phases and the computing effort has been recorded for 
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the complete formulation. In reference [5 ], the factors of 
matrices B* and B* * or their inverses and certain quantities 
such as C or Z have been assumed to be available in the 
preparatory phase and then comparison has been made for the 
•computing efficiency of the solution phase for the different 
compensation methods. Thus, in order to choose the best 
compensation method for a given application, a clear under- 
standing of compensation fundamentals and efficient implemen- 
tation is required. Skilful programming of sparse vector and 
matrix operations is also indispensable. 

The results of contingency evaluation for outage of 
line number 73 for I^EE 57 Bus system for Post, Pre and Mid- 
compensation based on hybrid approach are tabulated in Table 
4.5. These results obtained after one phase angle iteration 
and one voltage iteration on the base case load flow are quite 
compatible. 

4.4 CONCLUSION: 

Compensation methods can be applied advantageously 
over a wide range of applications, of which contingency 
evaluation is one such area. In this chapter, application 
of sparse vector methods has been examined in some of the 
. compensation techniques employed for contingency evaluation. 



Table 4 # 5: Results of Contingency Evaluation for IEEE 57 Bus System 


53 




0 

03 ^ 
*P • 


! 

t 

i 


















' | 




Q 

a 

a 

NO 

o 

o 

rH 

o 

o 

NO 

MT 

Q 

CM 

CM 

rH 

o 

CM 

Q 

n 





q 

o 

a 

rH 

lO 

o 

i — 1 

o 

o 

r- 

O 

O 

»H 

ON 


CO 


n 




O • 


a 

Q 

o 

ON 

00 

o 

o 

o 

a 

n 

00 

CD 

a 

rH 

CO 

n 

sD 

ON 

co 


c 

> Q- 

•H 

Q 

o 

lO 

rH 

NO 

a 

in 

n 

a 

ON 

n 

n 

«H 

n 

ON 

co 

r- 

o 

00 


o 

CO C 

> 


rH 

00 

CO 

l> 

00 

00 

Q 

00 

00 

r- 

rH 

00 

r- 

00 

rH 

rH 

rH 

a 


•H 

s — ' 

o 

o 

ON 

ON 

On 

ON 

ON 

O 

ON 

ON 

ON 

O 

ON 

ON 

ON 

o 

o 

o 

a 


-P 

D *H 


• 

• 

• 

• 

• 

• 

• 

• 

* 

© 

• 

© 

© 

© 

© 

* 

• 

« 

• 


03 

CO 

r-» 

CQ — 


rH 

rH 

o 

o 

Q 

Q 

a 

rH 

o 

o 

a 

H 

o 

a 

a 

rH 

•H 

rH 

rH 


CD 

a 























s 























O 

0 0 


o 









CM 

CM 

CM 






H 

a 


o 

■H 0 


tO 

ON 

NO 

rH 

CO 

rH 

NO 

a 

CD 

n 

l> 

VO 


n 

o 

ON 

00 

a 


1 

£T> 0 


o 

CM 

o 

CM 

CM 

NO 

rH 

ON 

H 

N0 

r- 

H 

rH 

ON 

rH 

H 

n 

n 

co 


T3 

C M 


o 

CM 

r- 

CO 

tO 

o 

CO 

a 

CD 

CM 

rH 

ON 

CD 

00 

O 

00 

n 

CM 

in 



CO O' 


Q 

if) 

lO 

CM 

sf 

ON 

in 

CO 

n 

O 

H 

00 

CO 

ON 

ON 

«H 

CM 

H 

n 


s 

0 

> 

Q 

C30 

D- 

CM 

CM 

ST 

00 

r- 

ON 


CM 

Nt 

CM 

ID- 

ON 


o 

n 

co 



0 T3 


O 

rH 

On 

CO 

in 

NO 

lO 


n 

• 

« 

* 

00 

CO 

rH 

00 


• 

© 



if) 


« 

* 

♦ 

© 

© 

• 

# 

♦ 

♦ 

rH 

CD 

o 

© 

• 

© 

• 

© 

H 

CO 



0 C 


o 

rH 

in 

o 

00 

oo 

0- 

'sT 

ON 

*H 

rH 

rH 

ON 

ON 

r- 

00 

n 

H 

H 



JC *h 
P k'-" 



l 

i 

I 

i 

1 

i 

I 

I 

1 

l 

1 

I 

1 

i 


i 

1 

I 



0 

CT> 

0 


o 




















Cf-P * 


Q 

o 

NO 

CO 

Q 

CO 

O 

o 

CO 

n 

Q 

CM 

M0 

O 

Q 

CM 

o 

o 


o 

1 » D 


o 

O 

a 

rH 

CO 

O 

in 

o 

a 

CO 

rH 

CD 

n 

M0 

CO 

00 

’sD 

if) 

rH 



O • 

> 

a 

O 

o 

ON 

CM 

o 

CO 

o 

a 

n 

rH 

O 

00 

00 


if) 

sD 

ON 



-P 

> a 

*H 

o 

a 

tn 

rH 

c- 

o 

i> 

n 

a 

CO 

t> 

n 

nO 

NO 

00 

CM 

ID- 

o 

Nt 


03 


"W 


rH 

00 

00 

o 

00 

00 

o 

00 

ON 


rH 

t> 

ID- 

00 

rH 

rH 

«H 

Q 


CO 

0 a 


o 

o 

ON 

ON 

ON 

ON 

ON 

o 

ON 

ON 

ON 

a 

ON 

ON 

ON 

D 

O 

o 

O 

* 

c 

n *h 


* 

• 

• 

• 

• 

• 

© 

• 

© 

© 

© 

© 

« 

• 

© 

© 

© 

© 

• * 


0 

01 
c; 

CX5 - — 


rH 

rH 

a 

Q 

a 

o 

o 

H 

o 

a 

o 

rH * 

Q 

o 

o 

rH 

rH 

rH 

»H 

DC 

O 

a 






















o 

i 

0 0 





















< 

0 

rH 0 




On 







CM 

CM 

CM 

nO 


n 

CD 

ON 

rH 

o 

o 

H 

O' o 

N 

O 

lO 

NO 

rH 

CO 

1 — 1 

NO 

a 

(D 

n 

O 

rH 

ON 

rH 

«H 

n 

00 

o 

05 

a* 

C H 


o 

CM 

CD 

CM 

CM 

NO 

rH 

ON 

rH 

nO 

r- 

H 

O 

00 

ON 

00 

n 

n 

CO 

A« 


< O' 

•H, 

o 

CM 

{> 

CO 

Ln 

t> 

CO 

o 

a 

CM 

rH 

O' 

CO 

ON 

ON 

«H 

CM 

CM 

n 

% 


0' 

«H 

o 

if) 

lO 

CM 

sr 

On 

in 

CO 

n 

t> 

«H 

00 

CM 

l> 

On 

I> 

Q 

«H 

in 

< 


0 T3 


o 

00 . 

tD- 

CM 

CM 

sr 

00 

o 

ON 


CM 


00 

co 

H 

00 


n 

co 



0 


o 

rH 

ON 

CO 

in 

nO 

n 

sr 

n 

.© ' 

© 

© 

© 

© 

• 

© 

© 

© 

• 

a 


fD c 


• 

• 

• 

© 

© 

© - 

« 

© 

© 

i — i 

o 

o 

ON 

ON 

r- 

CO 

n 

rH 

co 

i 


x: *h 


o 

rH 

in 

o 

00 

00 

o 

^ • 

ON 

rH 

rH 

*H 




H 

rH 

o5 

£p 


04'-' 



1 

i 

1 

I 

1 

i 

i 

l 

I 

1 

l 

I 

l 

1 

I 

i 

1 

1 

X 
























c 

0 






















o 

o> 























0 






















-P 

-P • 


p 

o 

o 

N0 

CO 

o 

CO 

o 

o 

CO 

n 

o 

CM 

NO 

CM 

CM 

CM 

a 

ID- 


03 

rH 0 


o 

o 

Q 

rH 

00 

O 

n 

o 

o 

CO 

rH 

o 

if) 

NO 

CO 



n 

r- 


03 

O • 


o 

o 

O 

ON 

CM 

O 

CO 

o 

o 

n 

rH 

o 

00 

00 


sD 


ON 



C 

> 0* 

*H 

o 

o 

in 

rH 

o 

CO 

i> 

n 

o 

00 

r- 

n 

nO 

NO 

00 

l> 

I> 

o 



0 


♦H 


rH 

CO 

CO 

r- 

On 

CO 

p 

00 

ON 

r- 

rH 

l> 


00 

rH 

«H 

rH 

o 


a 

0) C 


o 

o 

ON 

ON 

ON 

ON 

ON 

o 

ON 

ON 

ON 

o 

ON 

ON 

ON 

Q 

o 

Q 

o 


£ 

0 -H 


© 

• 

• 

© 

• 

© 

» 

■ • 

' • 

' # 

* 

• 

« 

© 

©. 

♦ ■ 

• 

© 

• 


o 

O 

CQ'— ' 


rH 


o 

o 

o 

O 

o 

rH 

o 

o 

o 

*H 

o 

o 

o 

rH 

rH 

«H 

rH 


i 

-p 























CO 

0 0 






















o 

rH u> 






















a* 

D> 0 





NO 






CM 

CM 

CM 






rH 

o 



B N 


o 

n 

ON 

1 — \ 

CO 

rH 

NO 

o 

o 

n 

i> 

NO 


in 

a 

ON 

00 

o 



<C O' 


8 

CM 

o 

CM 

CM 

NO 

•H 

ON 

rH 

VO 

i> 

rH 


ON 

H 

rH 

n 

n 

co 



o 

/^N 

CM 


CO 

n 

ID- 

CO 

o 

Q 

CM 

rH 

ON 

o 

00 

JN 

00 

n 

CM 

n 



0 ) -a 


Q 

lO 

lO 

CNi 

sr 

ON 

n 

CO 

n 

r- 

rH 

00 

CO 

ON 

ON 

rH 

CM 

rH 

n 

. 


V) 

'w' 

o 

00 

o 

CM 

CM 

sT 

00 

c^ 

ON 


CM 


CM 

ID- 

ON 

t"- 

a 

If) 

co 



o c 


o 

rH 

ON 

CO 

in 

N0 

n 

NnI* 

n 

© 

© 

« 

00 

CO 

rH 

00 


© 

© 



x: -h 


* 

• 

• 

• 

* 

© 

© 

© 

© 

* — i 

o 

o 

' © 

© 

© 

© 

• 

rH 

co 



lO*'—' 


o 

rH 

in 

r- 

00 

CO 

c- 

<sf 

ON 

H 

rH 

«H 

ON 

ON 

r- 

00 

n 

rH 

H 






1 

i 

I 

i 

t 

i 

i 

I 

1 

I 

1 

1 

1 

i 

r 

i 

I 


• 























»H 

O 



rH 

CM 

CO 

si- 

in 

NO 


00 

ON 

o 

rH 

CM 

CO 


n 

nO 

i> 

00 

0N * 

1 tn x 


V 










rH 

rH 

H 

rH 

rH 

rH 

rH 

rH 

«H 

-• 


« 

u 

-p 

c 

o 

o 

lO 

0 

H 

X> 

0 

H 




13.804270 1.012948 -13.804270 1.012840 -13.804270 1.020408 


lO 

rH 

O 

St 

St 


VO 

o 

CM 

ON 

o 

ON 

n 

in 

n 

CM 

o 

o 

CO 

St 

CM 

vO 

M0 

MO 

CO 

00 

MO 

St 

st 

rH 

St 

00 

MO 

ON 

CM 

MO 

On 

o 

St 

o 

CO 

i — j 

O 

<H 

st 

M0 

st 

rH 

st 

00 

rH 

CM 

CM 

00 

I> 

st 

CM 

r- 

00 

CM 

O 

CM 

CO 

co 

CO 

n 

CO 

\0 

CO 

CM 

rH 

l> 

o 

rH 

H 

\D 

MO 

CO 

00 

o 

rH 

n 

CM 

co 

«H 

M> 

o 

00 

On 

rH 

o 

O 

O 

o 

0 s 

On 

o 

o 

O 

ON 

ON 

ON 

ON 

O 

ON 

ON 

ON 

o 

# 

ft 

*■ 

• 

♦ 

« 

♦ 

. ft 

• 

♦ 

* 

• 

• 

♦ 

« 

• 

• 

• ft 

■ — 1 

*H 

rH 

o 

o 

o 

o 

rH 

•"H 

o 

o 

o 

o 

o 

o 

o 

O 

rH 



st 

ON 

CM 

MO 

*H 

CM 

ON 

in 

\0 

CM 

in 

O 

o 

st 

ON 

CM 

CM 

CM 

in 

in 

00 

«H 

rH 

00 

M0 

co 

rH 

O 

rH 

ON 

o 

CM 

00 

st 

M0 

ON 

ON 

vO 

00 

ON 

CO 

st 

n 

in 

n 

CM 

co 

St 

co 

vO 

o 

c- 

o 

00 

St 

ON 

H 

n 

in 

CM 

St 

On 

co 

O 

lO 

00 

o 

!> 

On 

ON 

I> 

St 

CM 

00 

O 

rH 

co 

CM 

vO 

CO 

st 

CM 

lO 

On 

o 

CM 

1> 

o 

CM 

CM 

00 

00 

ON 

CM 

i — f 

ON 

St 

st 

l> 


CO 

CO 

st 

rH 

00 

MO 

st. 

O 

ft. 

. ft 

ft 

• 

• 

*' 

• 

• 

• 

' * 

* 

e 

• 

♦ 

• 

ft 

* 

• 

CM 

CM 

CM 

co 

00 

CM 

rH 

o 

ON 

00 

ON 

00 

00 

Sf 

co 

CO 

CO 

CM 

rH 

«H 

rH 

H 

rH 

rH 

rH 

1 — { 


rH 

rH 

rH 

rH 

«H 

rH 

rH 

rH 

rH 

1 

l 

1 

1 

I 

l 

l 

I 

I 

s 

1 

I 

1 

1 

1 

J 

1 

l 



00 

St 

rH 

*H 

MO 


o 

in 

«H 

st 

00 

o 

00 

in 

in 

o 

CM 

ON 

rH 

CM 

MO 

f> 

ON 

co 

CO 

rH 

00 

st 

o 

00 

00 

M) 

in 

00 

St 

ON 

CO 

CM 

O 

ON 

st 

rH 

CM 

rH 

o 


o 

o 

St 

st 

St 

00 

c- 

M0 

CM 

O 

ON 

MO 

ON 

St 

MO 

00 

st 

CM 

CO 

in 

o 

CO 

CM 

• — { 

t> 

ON 

O 

O 

i> 

in 

l> 

O 

O 

rH 

co 

o 

o 

o 

n 

MO 

F- 

00 

o 

ON 

O 

O 

ON 

on 

ON 

O 

O 

O 

o 

o 

ON 

Ox 

0 s 

On 

ON 

Ov 

o 

• 

• 

• 

ft 

• 

ft 

ft 

« 

ft 

ft 

• 

ft 

ft 

ft 

•• 

# 

• 'ft 

ft 

o 

rH 

■«H 

o 

o 

o 

o 

rH 

rH 

o 

o 

o 

o 

o 

o 

o 

o 

rH 


St 

ON 

CM 

MO 

rH 

CM 

ON 

n 

M0 

CM 

n 

o 

O 

St 

Ox 

CM 

CM 

CM 

in 

Ln 

00 

rH 

«H 

00 

MO 

co 

i — | 

I> 

rH 

ON 

o 

CM 

00 

st 

MO 

O^ 

ON 

MO 

CO 

O 

CO 

st 

n 

n 

n 

CM 

co 

st 

co 

vO 

o 

!> 

o 

00 

St 


«H 

in 

in 

CM 

st 

ON 

co 

O 

n 

CO 

o 

t> 

ON 

On 

o 

st 

CM 

00 

o 

«H 

co 

CM 

MO 

CO 

st 

CM 

n 

ON 

o 

CM 

!> 

O 

CM 

CM 

00 

00 

ON 

CM 

rH 

Ov 

st 

st 

!> 

t> 

co 

co 

St 

rH 

00 

M0 

St 


ft 

ft 

ft 

# 

ft 

a 


ft 

• 

• 

- « 

• .. 

ft 

ft 

ft 

• 

* 

♦ 

CM 

CM 

CM 

00 

00 

CM 

rH 

o 

O' 

CO 

O' 

00 

00 

st 

CO 

CO 

co 

CM 

i — i , 

rH 

rH 

rH 

H 

rH 

rH 

rH 


« — ? 

rH 

rH 

* — 1 

rH 

«H 

rH 

rH 

rH 

1 

l 

1 

1 

1 

1 

I 

! 

i 

I 

I 

l 

I 

1 

1 

I 

i 

I 


00. CO 'O lO OJ H- CO iH O 00 VO ■ 00 O 

/Y"\ /T\ ■ ■■■ » /~v I «p 4« I . /V\ « — 4" U . — L * — W I 


CO 

in oo on 
st co 


H CM ^ H CO 
O (N 'O CM CO 


‘ » >4 \ / W Ov v J • » W V.M 

t" CM O O' r- O' ^ 00 in 04 00 tf) co •sf 

/-tv r~\ , r\ rv_ /-tv i ,-v-v /--> /— s . /-v . /-> 


O' O O l> iO I s O' O '— I 

O' O O O' O' O' O' O o 

1 . » • » . « • . . 


. __ _ „ _ -q- vo o vo 

■H OCO O 'O CM CM CM f- O' 

00 -H O 

C" 00 o 


00 O O O to 'O 


O'OO'O'0'OOO' O 


o 

rH 

iH 

O 

O 

o 

o 

rH 

' — f 

o 

o 

o 

O 

o 

o 

o 

o 

rH 

St 

ON 

CM 

M> 

rH 

CM 

ON 

in 

MO 

CM 

n 

o 

Q 

st 

ON 

CM 

CM 

CM 

in 

n 

00 

«H 

»H 

oo 

M) 

co 

rH 

t- 

rH 

ON 

o 

CM 

CO 

st 

M0 

O 

ON 

vO 

00 

ON 

co 

St 

n 

to 

n 

CM 

CO 

St 

co 

xO 

o 

i> 

o 

00 

St 

ON 

rH 

n 

n 

'M 

st 

ON 

co 

O 

n 

CO 

n 

t- 

ON 

On 

c- 

st 

CM 

co 

O 

rH 

00 

CM 

M) 

co 

st 

CM 

n 

ON 

o 

CM 

i> 

O 

CM 

CM 

00 

00 

ON 

■ CM 

rH 

ON 

St 

St 

cs 

t> 

v0 

co 

st 

H 

00 

NO 

st 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

.ft ' 

' 'ft' 

* 

ft 

* 

ft . 

• 

• 

CM 

CM 

CM 

00 

00 

CM 

•H 

o 

ON 

00 

ON 

CO 

CO 

St 

CO 

co 

co 

CM 

rH 

» — 1 

rH 

#H 

« — 1 

rH 

rH 

rH 


rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

1 

i 

I 

l 

r 

i 

I 

l 

I 

I 

1 

1 

J 

1 

1 

I 

I 

1 

rH 

CM 

co 

st 

n 

nO 

O 

CO 

ON 

o 

rH 

CM 

00 

st 

n 

M> 

l> 

00 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

co 

co 

CO 

CO 

co 

co 

CO 

CO 

co 


54 


Table 4.5: contd 
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Based on the cases studied, it is observed that a 
hybrid approach having a judicious mix of FF/FB 
approach for sparse vectors and Full Forward/Full Backward 
approach for non-sparse vectors^ may be more beneficial 
for network solutions^ 



CHAPTER 5 


CONCLUSION 

5.1 GENERAL: 

Notwithstanding, the relatively less exploitation 
of vector sparsity vis-a-vis the matrix sparsity, so far, 
there exists sufficient basis for gainfully employing sparse 
vector methods, to further enhance the efficiency of matrix 
solution methods. Sparse vector approach would be parti- 
cularly suited in power system applications where the inde- 
pendent vector B is sparse or only a subset of the solution 
vector X is required to be evaluated. Occurrence of sparse 
independent vectors is not uncommon in varied power system 
applications and thus sparse vector methods can lend them- 
selves to utilization in such areas as short-circuit analysis, 
optimal load flow, economic despatch, contingency analysis, 
compensation techniques etc. 

Normally in the solution of network equations, the 
Full Forward/Full Backward solution approach is resorted 
to. The sparse vector method based on the Fast Forward/ 
Fast Backward (FF/FB) approach is examined in Chapter 2, 
for solution of network equations. In several power system 
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applications, where the independent vector B is sparse or 
only a subset of the columns of vector X is to be determined, 
all nonzero elements are not required in the computation 
process for 'obtaining solution of linear algebraic equations. 
Factorisation paths can be determined using algorithms 
described in Chapter 2 for identifying the relevant nonzero 
entries in the sparse A matrix,for the solution of network 
equations. 

A comparison of Full Forward/Full Backward approach 
with FF/FB approach suggests that Full Forward/Full Backward 
method may be preferred when the vector B is full or contains 
large number of nonzero elements. This is because when 
vector B is predominantly non-sparse, the path finding algorithm 
becomes time consuming and consequently inefficient. However, 
the FF/FB approach is faster and hence suitable, when there 
are very few nonzero entries in vector B or when vector B is 
a singleton. The computational efficiency of this technique 
is further enhanced with increasing sparsity of system 
matrix A. 

Having turned out to be effective in the solution of 
linear algebraic equations, when the independent vector is 
highly sparse, the viability of sparse vector approach in 
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power system applications viz. contingency evaluation has 
been explored. Based on case studies reported in Chapter 3, 
it may be commented that neither the FF/FB method (i.e. 
the sparse vector approach) nor the Full Forward/Full 
Backward approach is beneficial by itself. A hybrid approach 
integrating FF/Fb method with Full Forward/Full Backward 
method, exploiting their individual merits has proved to be 
more effective, ■ 

Compensation techniques are well suited for contin- 
gency evaluation. In Chapter 4 , investigation of three 
compensation techniques namely Post, Pre and Mid-compensation 
in contingency evaluation, with a view to determining the 
feasibility of sparse vector approach suggests that the sparse 
vector approach is capable of improving the computational 
efficiency of these compensation methods. 

Case studies on various sample power systems seem to corro- 
borate the observation that the hybrid approach would 

be more competitive, than either the FF/FB or the Full Forward/ 
Full Backward method. Moreover, though the relative effecti- 
veness of these compensation techniques in contingency evalua- 
tion has also been analysed, a detailed study of this aspect 
would be necessary before coming to any definite conclusions , 
in this regard. 
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5.2 SCOPE FOR FURTHER WORK: 

While the application of sparse vector approach to systems upto 
100 buses has been examined, it is important to examine the 
effectiveness of this approach over larger power system 
networks. 

The effectiveness of the sparse vector approach also 
needs to be investigated when sparsity oriented storage 
schemes are employed for sparse matrix A and optimal ordering 
is resorted to. 

It would also be relevant to study the effectiveness 
of the hybrid approach for contingency evaluation, in greater 
depth over large systems, considering multiple simultaneous 
outages. 
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